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The boundary-induced phase transitions in an asymmetric simple exclusion process with inter- 
particle repulsion and bulk non-conservation are analyzed through the fixed points of the boundary 
layers. This system is known to have phases in which particle density profiles have different kinds 
of shocks. We show how this boundary-layer fixed-point method allows us to gain physical insights 
on the nature of the phases and also to obtain several quantitative results on the density profiles 
especially on the nature of the boundary-layers and shocks. 



I. INTRODUCTION 

Asymmetric simple exclusion process (ASEP) is a non- 
equilibrium process in which particles hop to the neigh- 
boring site in a specific direction following certain con- 
ditions In the simplest model, particles obey the 
mutual exclusion condition due to which a site cannot 
be occupied by more than one particle and a hop to the 
neighboring site is possible provided this target site is 
empty. The particles after being injected at one end of 
the lattice at a rate a hop across the lattice and reach 
the other end where they are withdrawn at a rate 1 — 7. 
There exist other models where particles can have attrac- 
tive or repulsive interaction in addition to the exclusion 
interaction @,0|. 

All these systems exhibit interesting boundary-induced 
phase transitions for which the tuning parameters are the 
boundary rates, a and 7 [HHf. In different phases, the av- 
erage particle density has distinct constant values across 
the bulk of the lattice. These particle density profiles 
in different phases may also differ due to different loca- 
tion and nature of the boundary-layers. More features, 
such as coexistence of high and low-density regimes are 
seen in systems where particle number is not conserved 
in the bulk due to attachment and detachment of par- 
ticles to and from the lattice [1, • in this coexistence 
phase (also known as a shock phase), the particle density 
profile has a jump discontinuity (shock) in the interior 
of the lattice from a low to a high density value. The 
slope and the number of such shocks in a density profile 
are related to the nature of the inter-particle interactions 
which determine the fundamental current density rela- 
tion [1, Although this relation predicts the kind of 
shocks that can be seen Q in the density profile, a sys- 
tematic characterization of the phase-transition and the 
phase diagram in the a — 7 plane requires a detailed anal- 
ysis of the relevant equations describing the dynamics in 
the steady-state @, OH EH • Drawing analogies from the 
equilibrium phase transitions, first order, critical 0, @] 
and tricritical [ll| kind of phase transitions have been 
observed so far in various ASEP models. 

While characterizing the phase transitions, it is useful 
to study the variation of the height or the width of the 
boundary-layers as a and 7 are changed. In a way, these 



boundary- layers play important roles in deciding the "or- 
der parameter" like quantities in these non-equilibrium 
phase transitions. Owing to their importance in describ- 
ing phase transitions, the boundary-layers for several in- 
teracting and non-interacting models have been studied 
using the techniques of boundary- layer analysis (l2j . It is 
found that the system usually enters into a shock phase 
from a non-shock phase due to the deconfinement of the 
boundary-layer from the boundary. This deconfinement 
can be described by a nontrivial scaling exponent asso- 
ciated with the width of the boundary- layer Q. For the 
pure exclusion case, except for a critical point, the transi- 
tion from a non-shock to a shock phase is first order in na- 
ture since a shock of finite height is formed on the phase 
boundary. The height of the shock on the phase bound- 
ary reduces as one approaches the critical point along the 
phase boundary. At the critical point, the shock height is 
zero on the phase boundary and it increases continuously 
as one proceeds away from the phase-boundary further 
into the shock phase. In order to visualize these features, 
it is beneficial to obtain the full solution for the density 
profile along with its boundary-layer. Boundary layer 
analysis is useful for this purpose since it allows us to 
generate a uniform approximation for solving the steady- 
state particle density equation across the entire lattice. 
This steady-state equation can be obtained from the large 
time- and length-scale limit (hydrodynamic limit) of the 
statistically averaged master equation that describes the 
particle dynamics in the discrete form. For the simple 
exclusion case, it is possible to obtain an analytical solu- 
tion of the steady-state hydrodynamic equation for the 
entire density profile. This, however, may not possible 
for more complex ASEPs. 

A fixed-point analysis of the hydrodynamic equation 
turns [l3| out to be general and useful, since it docs 
not involve an explicit solution of the steady-state hy- 
drodynamic equation. In particle conserving models, a 
boundary-layer saturates to the constant bulk density 
profile asymptotically. As a consequence of this, it is ex- 
pected that the fixed-points of the boundary-layer equa- 
tion match with the bulk density values. In other words, 
a boundary-layer, which is a solution of the boundary- 
layer equation, is a part of the flow trajectory of the equa- 
tion flowing to the appropriate fixed-point on the phase 
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plane. Thus, in order to find out the values of the bulk- 
densities in different phases, it is sufficient to determine 
the physically acceptable fixed-points of the boundary- 
layer equation. As a result, the number of possible bulk 
phases is given by the number of these fixed-points. Ap- 
plying this method to a specific particle conserving two- 
species ASEP [l3| , it is found that this system has three 
distinct bulk phases corresponding to three fixed-points 
of the boundary-layer equations. In addition, it is possi- 
ble to predict the nature of phase transitions, locations 
of the boundary-layers etc. for this system. All these 
predictions match well with the results from numerical 
simulations [Til ]. 

In a particle non-conserving case, the density is not 
constant in the bulk, and therefore, the fixed-points of 
the boundary-layer do not provide the full profile since 
the details of the bulk dynamics is not considered in 
this approach. However, it is still useful to obtain the 
fixed-points of the boundary-layer equations along with 
their stability properties in order to predict the possible 
shapes of the density profiles under different boundary 
conditions. In the present paper, we consider a particle 
number non-conserving model where particles interact re- 
pulsively. Our aim is to extend the fixed-point analysis 
to a system with non-constant bulk density. We show 
how this analysis helps us predict possible shapes of the 
density profiles under different boundary conditions and 
also understand the properties of different kinds of shocks 
present in the density profile. This particular model is 
chosen because of certain nontrivial shapes of the density 
profiles with different kinds of shocks. 

The plan of the paper is as follows. In the following 
section, we describe the model. This section also con- 
tains brief discussions on the hydrodynamic approach, 
boundary-layer analysis and some of the known results. 
In section III, we present the phase-plane analysis of the 
boundary-layer equations for the present model. There 
are separate subsections on the boundary-layer equation, 
its fixed-points, stability analysis of the fixed-points and 
possible shapes of shocks. Section IV presents the predic- 
tions of the possible shapes of the density profile under 
different boundary conditions. In section V, we mention 
a few general rules for predicting the shapes of the density 
profiles and some special features related to the shocks of 
this model. We end the paper with a summary in section 
VI. 



II. MODEL 

A. Discrete description 

The asymmetric simple exclusion process that we con- 
sider here consists of a one-dimensional lattice of N sites 
with lattice spacing, a. Particles are injected at i = 1 
with rate a and withdrawn at i = N at a rate 1 — 7. Par- 
ticles, obeying mutual exclusion, hop to the right with 
rates that depend on the occupancy of the neighboring 



site as 



1100 1010 at a rate 1 + e, 

0101 0011 at a rate 1 - e, 

0100 0010 at rate 1 

1101 1010 at rate 1. 



(1) 
(2) 
(3) 
(4) 



Here, < e < 1 and 1 (0) represents an occupied (un- 
occupied) site. For e 7^ 0, there is an effective repulsion 
between the particles @, H, Q . In addition, the number 
of particles is not conserved due to particle detachment, 
1 — > 0, at a rate uJd and attachment, — > 1, at a rate uj a 
at any site on the lattice. Particle attachment and de- 
tachment are equilibrium like processes that do not give 
rise to any particle current. 



B. Hydrodynamic Approach and a brief 
description of the boundary-layer analysis 

The hydrodynamic approach is based on the lattice 
continuity equation which equates the time evolution of 
the particle occupancy at a given site with the difference 
of currents across its two neighboring bonds. In the con- 
tinuum description, the continuous time and space vari- 
ables are t and x with the latter replacing, for example, 
the ith site as i — > x = ia. Upon doing a Taylor ex- 
pansion of the statistically averaged continuum version 
of the lattice continuity equation in small a, one has the 
following hydrodynamic equation 



dp + dJ + S -0 

-57 + T\ H o — U, 

at ox 



(5) 



for the averaged particle density p(x,t). This equation 
has already been supplemented with the particle non- 
conserving parts 



5*0 = -&{pl - p), 



(6) 



where pl ~ 



and f2 = {u a + LOd)N. The current, 



J(p), consists of a bulk current j{p) and a diffusive cur- 
rent proportional to as 



J 



9P _L ■( \ 



(7) 



Here, eo is a small parameter proportional to a. The dif- 
fusive current part arises naturally as one retains terms 
up to 0(a 2 ) in the Taylor expansion. In order to deter- 
mine the particle density, p(x), in the steady-state (§f = 
0), one has to solve the differential equation with appro- 
priate boundary conditions. We consider the lattice-ends 
to be attached to the particle-reservoirs which maintain 
constant densities p{x = 0) = a and p(x = 1) = 7. 
The diffusive current part is crucial here since due to its 
presence, the hydrodynamic equation becomes a second 
order differential equation and as a result we can obtain a 
smooth solution satisfying both the boundary conditions. 
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e = is the usual ASEP with only the exclusion 
interaction. In this case, the current density relation, 
j(p) = p{\ — p) is an exact one. The symmetric shape 
of the current about its maximum at p = 1/2 is a con- 
sequence of its invariance under particle-hole exchange 
p — > 1 — p. It is well understood that the phase dia- 
gram has low-density (p < 1/2 at the bulk), high-density 
(p > 1/2 at the bulk) and maximum current (p = 1/2 at 
the bulk) phases Q. The particle- hole symmetry is re- 
tained in e 7^ models although the current changes non- 
trivially. At e = 1, i.e. for the extreme repulsion case, 
hops such as 0101 — > 0011 are forbidden. The current, 
therefore, vanishes exactly at the half-filling (p = 1/2) 
with the maximum current appearing symmetrically for 
densities on the two sides of p = 1/2. The exact form 
of the current as a function of p for arbitrary e can be 
found using a transfer matrix approach [U and it evolves 
from a single to a symmetric double peak structure as e 
grows beyond e,/ w .8. A simple, analytically tractable 
form of the current with double peaks can be obtained 
by doing a double expansion of the exact current about 
e = e,/ and p = 1/2 [11]. This leads to a quadratic form 
for the current 

j(p) = (2r + u)/16 - T -(p 1/2) 2 - u(p 1/2) 4 , (8) 

where the constant term is chosen in such a way that 
j(p) = for p = 0, or 1. We recover the non-interacting 
limit, j(p) = p(l — p) for r = 2 and u = 0. The double 
peak shape appears for r < 0. In the entire analysis 
below, we consider r to be a small negative parameter 
and u > 0. 

For the boundary-layer analysis, it is important to con- 
sider the bulk part and the narrow boundary-layer or 
the shock regions of the density profile separately. These 
boundary-layers or shocks are formed over a narrow re- 
gion of width O(fo) an d they merge to the bulk den- 
sity in the appropriate asymptotic limit. In order to 
study the boundary-layer and its asymptotic approach to 
the bulk, one can rescale the position variable in §5§ as 
x = (x — xq) / eo, where Xq is the location of the center of 
the boundary-layer. Hence, for a boundary-layer satisfy- 
ing the boundary condition at x = 1 , we have xq ~ 1 . For 
small eo, the boundary-layer approaches the bulk density 
in the x — > — oo limit and satisfies the boundary condition 
at x = 0. In terms of x, the steady-state hydrodynamic 
equation is 



(Pp 

dx 2 



dx 



e Q S Q = 0. 



(9) 



Since eo is a small parameter, the effect of the particle 
non-conserving term, So, on the boundary-layer is neg- 
ligible. As a result the total current, J = j(p) — || is 
constant across the boundary-layer. A shock, therefore, 
can be represented by a horizontal line connecting two 
densities in the j — p plane as shown in figure ([T]). For an 
upward shock > 0), this line lies below the j(p) curve 
and the reverse is true for a downward shock < 0). 




FIG. 1: Current j is plotted as a function of p. Lines 'a' and 
'b' represent an upward and a downward shock respectively. 



As a result, while for r > 0, only upward shocks are 
possible, for r < 0, there can be density profiles with a 
downward shock and double shocks. Double shocks can 
be represented by two horizontal lines on the j — p plane 
below the two peaks in j(p). 

To zeroth order in eo, the final boundary-layer equation 

is 



(Pp 
dx 2 



dx 



0. 



(10) 



In the boundary-layer language, the solution of this equa- 
tion is known as the inner solution. To obtain the bulk 
part of the density profile, one can ignore the diffusive 
current part in J for small eo- The steady-state equation 
that gives the bulk part of the density profile is 



dx 



Sn=0. 



(11) 



The solution of this equation for the bulk part of the pro- 
file is known as the outer solution. These inner and outer 
solutions contain several integration constants which are 
fixed by the boundary conditions and other matching 
conditions of the boundary-layer and the bulk under var- 
ious limits. Since the slope of the outer solution is ob- 
tained from 



dp 
dx 



So/f, 
dp 



(12) 



for a given p, the slope depends crucially on the signs of 
(pL — p) and For all the analysis below, we consider 
Pl to be large and a and 7 to be much smaller than p^. 



C. Known results 

Double peak structure of the current-density relation 
leads to two maximum current and one minimum cur- 
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rent phases in the phase diagram of the particle con- 
serving repulsion model Q. In the maximum and min- 
imum current phases, the bulk density values are those 
at which the current attains its maximum and minimum 
values respectively. With these new phases, the phase 
diagram for this model becomes more complex than its 
non- interacting counterpart. 

Combining the techniques of boundary-layer analysis 
and the results from numerical solutions, the phase dia- 
gram has been obtained for the particle non-conserving 
repulsion model 1_1[ . The phase diagram has a lot of in- 
teresting features including a tricritical point at r = 0. In 
the a — 7 — r phase diagram, this is a special point where 
two critical lines meet. It has been found that three dif- 
ferent phase diagrams are possible for r > 0, r = and 
r < 0. For r > 0, the current-density plot is symmetric 
around p = 1/2 with a maximum at p = 1/2. The nature 
of the phase diagram is qualitatively similar to the mu- 
tually exclusive case with one single critical point. For 
r < 0, with a double peak structure of the current-density 
plot, the phase diagram is more complex with more than 
one critical point and three different shock phases with 
the density profile having a single upward shock, double 
upward shocks and one upward and one downward shock 

m. 



The low-density peak can give rise to a low-density up- 
ward shock (p(x) < 1/2 in the shock part) in the density 
profile. A single shock of this kind can be represented 
by a horizontal line in the j — p plane below the low- 
density peak. The critical point corresponds to a situa- 
tion where the horizontal line reaches the peak position 
implying a shock of zero height. The second distinct 
critical point that involves both the peaks of the current- 
density plot is not symmetrically related to this. The 
density profile, here, has two upward shocks, in which 
one is a low-density shock and one is a high-density shock 
with p > 1/2. The low-density shock, in this case, has 
the maximum height with its high-density end saturating 
to p = 1/2. The high-density shock which is due to the 
high-density peak of the current-density plot can be of 
varying height. The critical point corresponds to the spe- 
cial point where this high-density shock has zero height. 
In addition to these regions, there are regions in the phase 
diagram, where density profiles with a downward shock 
or a single, symmetric upward shock arc found. 



III. PHASE-PLANE ANALYSIS OF THE 
BOUNDARY-LAYER EQUATIONS 

In the following subsections, we determine the fixed- 
points of the boundary-layer equation and their stability 
properties. These fixed-points are the special points to 
which the boundary-layer solution saturates in the ap- 
propriate limit. The knowledge about the fixed-points 
and their stabilities can, therefore, be used to our advan- 
tage to find out, for example, the bulk densities to which 
a shock or a boundary-layer saturates at its two edges. 



A. Boundary-layer equation 

Substituting the expression for j(p) as given in ([5} and 
integrating the boundary-layer equation, (fTOj) once, we 
have 



dpi 
dx 



"Pi i r 2 i u 4 n 

— + 4P1 + oPl = C °- 



(13) 



Here p\ = 2p—l and Co is the integration constant. The 
saturation of the boundary-layer to the bulk density, pu,, 
is ensured by choosing the integration constant as 



r 2 u 4 



(14) 



As per equation ((8j), Co is related to the excess current 
(positive, negative or zero) measured from p = 1/2 (half- 
filled case). The entire analysis in the following is done 
in terms of p\ for which the boundary conditions are 
P\{x = 0) = ai = 2a — 1 and p\{x = 1) = 71 = 27 — 1. 



B. fixed-points 

Co can be plotted for various pw varying from —1 to 
1. For r < 0, Co has a symmetric double well structure 
around pu = (see figure dH). The fixed-points, p*, of 
equation (|13|) , are the solutions of the algebraic equation 

(15) 



^ + ^-C„ = 0. 



In general, there are four possible solutions for the fixed- 
point as 



pi = ±[ ] 



±yfr 



8C 



0^11/2 



(16) 



In view of the symmetry of the j — p diagram, it is 
natural to expect the two critical points to be related 
through this symmetry. Previous work, however, shows 
that the shapes of the density profiles are not related 
through this symmetry near these two special points. Un- 
like the low-density shock, the high-density shock in the 
density profile is always accompanied by a low-density 
shock of maximum height. The following analysis clearly 
reveals the reasons behind such asymmetries. 



The value of Co depends on pu,, the bulk density to which 
the boundary-layer solution saturates. As a consequence, 
for a given Co, the corresponding pu is always a fixed- 
point. For the same Co, there are, however, other fixed- 
points which are determined from (|TB| . Hence, from the 
information about one saturation density pib, the other 
saturation density of the shock can always be determined. 
The approach to various fixed-points has to be, of course, 
consistent with their stability properties. These stability 
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FIG. 2: Co plotted as a function of pit with r 
u = 2.2 



.2 and 



properties of various fixed-points are discussed in the fol- 
lowing subsection. 

If Cq is positive, there can be only two real fixed-points 
of opposite signs. The positive and negative fixed-points 
denoted respectively as pl± are 



Pi± 



| r | +Vr 2 + 8C u in 



(17) 



If Co < 0, there are four fixed-points for pu,. In all these 
cases, the fixed-points are symmetrically located on the 
either side of the origin. The positive p\ fixed-points are 



Pl,2+ = I 



r | ±Vr 2 -8 | Co | u n 



/2 



and the negative p\ fixed-points are 



Pl.2- = 



±^ 2 -8|C 



ll/2 



(18) 



(19) 



Here, the subscripts 1 and 2 correspond to the + and 
— signs inside the square bracket respectively. It is im- 
portant to notice that for Co < 0, all the fixed-points 
become imaginary when | Co |> As Co approaches 
this lowest negative value, the pair of fixed-points on the 
positive and negative sides approach each other and they 

merge at Co = — At this special value, the fixed- 
points are p*^ = ±(-^) 1//2 . For Co = 0, there are three 
fixed-points, pi = and p*^ = i(^) 1 / 2 . 

Numerical values of the fixed-points for some special 
values of Co with r = — .2 and u = 2.2 are mentioned 
below. For Co = 0, the nonzero fixed-points are p*^ = 
0, ±.426. For these values of r and u, no real fixed-points 

2 

arc present if Co < — = —.00227. At this special value 
of Co, the two fixed-points are = ±.30151. 



C. Stability analysis of the fixed-points 



For Co > 0, a linearization of equation (|13j) around 
the fixed-points with pi = p\+ Spi leads to the following 
stability equation 



dSpi 
d.f 



-yV \ 2 +8uC 

o Pl°Pl- 



(20) 



This implies that the fixed-points p\ + and p\_ arc, re- 
spectively, stable and unstable. 

Similarly, for Co < 0, the general stability equation is 



d5pi 
d.f 



Pi 



(| r | - upf). 



(21) 



The flow around the fixed-points can be obtained by sub- 
stituting the explicit expressions of the fixed-points. Fig- 
ure © shows the stability properties of various fixed- 
points for Co > and Co < 0. 



7- 



c >o 



< — l — <- 





1+ 



2+ 



1 + 



FIG. 3: flow behavior of the fixed-points 



The stability property of the p* = fixed-point for 
Co = and for the pair of fixed-points for Co = — 
cannot be determined from the linear analysis. How- 
ever, the flow around the fixed-points can be predicted 
from the continuity of the flow behavior as Co approaches 
these special values. Fixed-points, their stability proper- 
ties and how the fixed-points change with Co are com- 
binedly shown in figure 



D. Shocks for various Co 

Since Co can be expressed purely in terms of J, with 
its constant parts subtracted, it remains constant across 
a shock or a boundary-layer. In principle, using equa- 
tion ([T4| . one can obtain the value of Co along the con- 
tinuously varying parts of the density profile. Hence, 
as we move along a density profile having bulk shocks, 
Co changes as per equation (fT4| along the outer solu- 
tion parts of the profile with intermediate constant values 
across the shock or inner solution regions. The value of 
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C () 

FIG. 4: Fixed-points are plotted for different values of Co 
with r = —.2 and u = 2.2. p*± and p%± are different fixed- 
points as mentioned in the text. The vertical solid lines with 
arrows show the flow behavior of the fixed-points. 



Co in the shock region is fixed by one of the bulk density 
values to which the shock saturates. 

With the information on the possible values of Co in 
the entire range of pu and, hence, the knowledge about 
the corresponding fixed-points, it is possible to list the 
kind of shocks that can be observed. 

Let us assume that the shocks or the boundary-layers 
approach the bulk densities p\ r or pu as x — > ±oo respec- 
tively. Since pu and p\ r are various fixed-points of the 
inner equation, the approach to these fixed-points has to 
be consistent with the flow properties. A shock is called 
an upward shock if pu < p\ r . The reverse i.e. pu > p\ r 
is true for a downward shock. 

(i) Co > 0: In this case, there are two fixed-points 
p* + and p\_ symmetrically located around p\ = with 
p*_ being an unstable fixed-point. Thus if a shock is 
formed with Co > 0, it should be an upward shock which 
approaches the fixed-points p\ r = p* + and pu = p*_ as 
x — > oo and — oo respectively. The shock height, in this 
case, is p* + — p\_. 

(ii) Co < 0: In this case, four fixed-points lead to dif- 
ferent kinds of shocks. 

(a) It is possible to see a downward shock with pu = 
pi, + and p\ r = p\_. The downward shock is thus sym- 
metric around p\ = 0. The flow in figure ([3]) shows that a 
downward shock cannot involve other fixed-points since 
that would not be consistent with the stability criteria of 
the fixed-points. 

(b) There can be small upward shocks which lie en- 



tirely in the range p\ > 0. We have already referred 
these shocks as high-density shocks. In terms of the fixed- 
points, the left and right saturation densities of the shock 
are pu = P2 + and p\ r = p* + , respectively. 

(c) The third possibility is that of an upward shock 
entirely in p\ < range. Such a shock has been referred 
as a low-density shock. For this shock, pu = p*_ and 

Plr = P2-- 

(iii) Co = 0: There can be an upward shock with p\ r = 
and pu = p\_. There can also be an upward shock 
connecting the densities pu = and p\ r = p* 1+ . These 
two shocks together appear as a large shock, symmetric 
around pi = 0. 

Alternatively, different kinds of shocks can tell us the 
range of values for Co . 



IV. PREDICTIONS ABOUT THE SHAPES OF 
THE DENSITY PROFILES 

Based on figure ([4]), we attempt to predict possible 
shapes of the density profiles for given boundary condi- 
tions oi\ and 71 . We consider only a few pairs of boundary 
conditions and based on this, we make certain general 
predictions in the next section. The basic strategy for 
drawing the density profile is as follows. We first need to 
mark a\ and 71 on the p\ axis of Co — p\ plane. Starting 
with cither of the boundary conditions, we change pi, 
along the curve in figure (|4]) , in a way that we reach the 
other boundary condition in the end of our move. While 
doing so, we may allow a discontinuous variation of p\ 
along a vertical constant-Co line, provided it does not 
violate the flow property. Such a discontinuous change 
in pi appears in the form of a shock or a boundary-layer 
in the density profile. The dashed, vertical lines in figure 
([5^,), for example, are the constant-Co lines along which 
the density may change. Such a dashcd-linc, therefore, 
represents a boundary-layer or a shock in the density pro- 
file. Two densities at which a shock or a boundary-layer 
saturates, are those at which a particular, constant-Co 
line, representing a shock or a boundary-layer, intersects 
the curves. This method, however, sometimes leaves us 
with different options for the density profile. All these 
possibilities are shown on Co — pi plane for each pair of 
boundary conditions individually. 



A. Density profiles with only boundary-layers 

Suppose we consider a situation where ati, 71 > with 
a\ < pl^ and 71 > p*^. There can be a possibility where 
the density profile has a particle-depleted boundary-layer 
|a; = o> 0) at x — satisfying the boundary condi- 
tions pi(x = 0) = oi\. This boundary- layer can be repre- 
sented by a vertical line similar to (a) in figure ((5ji). This 
is consistent with the flow property that suggests the ap- 
proach of the boundary-layer to the fixed-point p\ + as 
x — > 00. On the other hand, in the x — > —00 limit, 
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FIG. 5: (a) Possible variations of the density, as one moves 
along the density profile from x — end, are shown on the 
Co — Pi plane. Bold dashed lines or curves with open arrows 
show the variation of the density. Open arrows point in the 
direction of increasing x. Dotted lines mark the boundary 
conditions ot\ and 71. Solid lines with arrows are the flow 
trajectories of the fixed-points, (b) Numerical solutions for 
the density pi for various ai with 71 = 0.34. The inset gives a 
zoomed view of the particle-depleted boundary-layer at x = 0. 
No boundary-layer is formed near x = 1. 



which corresponds to the unphysical negative x region, 
the boundary-layer saturates to the unstable fixed-point 
pl_. After the boundary-layer, the density may decrease 
continuously along (c) on the p\ + branch and satisfies the 
boundary condition, pi(x = 1) =71. There can be an- 
other possibility where the particle-depleted boundary- 
layer at x = is represented by a vertical line similar 
to (b) joining the fixed-points p\ + and p* + . The bound- 
ary condition at x = 1 is again satisfied by a decreasing 
density part parallel to (c). For this to be possible the 
condition 71 < p*Q is required. These two possibilities 
are distinct due to distinctly different values of Cq . This 
shows the crucial role played by Co in deciding the den- 
sity profile. Numerical solutions of the full steady-state 
hydrodynamic equation presented in figure ^jp) show the 
boundary-layers saturating to a bulk density pu ~ .48. 
This implies that the boundary-layers are indeed repre- 
sented by (a) type vertical lines. 

(2) Next, we consider 01,71 > with a% > p*^ and 
71 < Pj^j. In this case too, the density profile can satisfy 
the boundary condition at x — through a boundary- 
layer that can be represented by a line similar to (a) in 
figure ([Sk). This would be a particle depleted boundary- 



layer at x = 0. In order to satisfy the other boundary 
condition, the density should decrease till along path 
(b) on the p\ + branch and then satisfy the right bound- 
ary condition through a particle-depleted boundary-layer 
along (c) in figure (|6^i). The boundary condition at x = 
can also be satisfied by vertical lines coming from above 
the p\ + branch leading to particle-rich boundary-layers 
U=o< 0) at x = 0. These lines could be (a') type 
lines in figure (|6jt) satisfying the boundary condition a[ . 
Both particle-depleted and particle-rich boundary-layers 
are present in the density profiles of figure ((BJd) obtained 
by solving the full steady-state hydrodynamic equation 
numerically. 



B. Density profiles with upward shocks 

With 71 > p*+ and ol\ large negative, possible shapes 
of the density profile can be of the following kinds. We 
start with p\(x = 0) = a\. p\ decreases continuously 
along p\_ branch along the dashed line (a) in figure ((7Ji). 
After this part, an upward shock, represented by a verti- 
cal line similar to either (b) or (d) appears. In the latter 
case, the dashed line (a) should be extended further till 
it reaches the low-density end of (d) . If the shock is rep- 
resented by (b), it is a large shock, symmetric around 
pi = 0. In the second case, the shock is a low-density 
shock with the saturation densities being p\ r = and 
P11 = P10 ■ If the density approaches the p* + branch af- 
ter the large shock, the boundary condition at x = 1 can 
be satisfied after that by a decrease in density along (c) 
on this branch. If the shock is of (d) kind, the density 
has to change further to satisfy the right boundary con- 
dition. Upon reaching pi = value, the density may 
change along p\_ or the pZ, + branch. The flow around 
p\ = 0, however, suggests that the density variation only 
along branch (path (c) in the figure ([7])) is possi- 
ble. The continuously increasing part along (e) is then 
followed by another upward shock, given by line (f ) , tak- 
ing the density to p* 1+ branch. The boundary condition 
is then satisfied by a continuously decreasing part along 
a (c) type line. Numerical solutions in figure (JTJj) are 
consistent with these predictions. 



C. Density profiles with downward shocks 

We next consider a case p\~ n < 71 < with ol\ in- 
creasing from large negative values. Here wc specifically 
mention how the density profile changes as ct\ is changed 
keeping 71 fixed. In the process, we observe how a den- 
sity profile with a downward shock appears. Let us as- 
sume that our starting ot\ lies somewhere on the p\_ 
branch, pi increases from p\{x = 0) = a.\ along p*„ 
branch till it reaches the boundary x = 1. This con- 
tinuously increasing part is represented by (la) in fig- 
ure (JHta). The density, then, satisfies the right boundary 
condition through a boundary-layer which can be, for 
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FIG. 6: (a) Same captions as figure J5K) except that here we explicitly show possible density variations for two different left 
boundary conditions specified by ot\ and ai. (b) Numerical solutions for pi (a;) for various ai with 71 = 0.23. The main figure 
provides a zoomed view of the boundary-layers near x = 0. The entire density profile over the entire lattice is shown in the 
inset, (c) A zoomed view of the same density profiles near the x = 1. This shows the particle depleted boundary-layers near 
x = 1. 



example, represented by a vertical line like (lb), in fig- 
ure (JSJi). This line takes the solution from the unstable 
fixed-point p\_ to the stable fixed-point p\ + . Since the 
vertical line (lb) passes through 71 before reaching the 
p* + branch, the boundary-layer satisfies the boundary 
condition before saturating to the positive fixed-point, 
p* 1+ . The boundary-layer at x = 1 is, therefore, a part of 
this vertical, constant-Co line. 

As «i is increased slightly, the route of the density 
along p*_ branch remains the same but this time the 
density reaches a higher value than the previous ot\ case 
before increasing sharply as a boundary-layer satisfying 
the boundary condition at x = 1. As a 1 is increased fur- 
ther, for a given ot\, the continuously increasing part of 
the profile reaches the low-density end of the line (2b). 
After this, there is a shock in the density profile of (2b) 
kind. This is a low-density shock that takes the density 
to p\ = 0. If this jump is near the boundary, this shock 
becomes actually a boundary-layer that can help the den- 



sity satisfy the boundary condition at x = 1. However, 
if this discontinuity is in the bulk, it is an upward low- 
density shock. In case of a shock in the bulk, the density 
increases further along pZ, + branch (path (2c) in figure 
([5])). The boundary condition, however, demands a de- 
crease in pi . This is possible through a downward vertical 
line (similar to path (2d)) and then a continuously de- 
creasing part (2e) along p\_ branch. The path (2d) is a 
downward shock that is seen in figure (|5Jd) and (|5]c). In 
case the density varies along lines (3c) and (3d), we have 
a downward boundary-layer near x = 1. These possibili- 
ties are expected if a\ is increased further from its value 
that leads to (2c) and (2d) type variations. As before, the 
principle is that if the (3d) type vertical line intersects 71 
line before reaching the p* 2 _ branch, (3d) type line repre- 
sents a boundary-layer satisfying the boundary condition 
at x = 1. If the reverse happens, this downward vertical 
line represents a downward shock at the bulk which needs 
to be followed by a continuously decreasing density part 



_0.6 ' ' ' ' ' ' ' ' ' ' 
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FIG. 7: (a) This figure has the same caption as figure Q. 
(b) Plot of the density profiles pi (x) for various large negative 
values of aj with 71 = .32. No boundary-layer is formed at 
x = or x = 1. A few density profiles towards right have 
double shocks. 



along the p^_ branch. This is a general principle which 
can be applied to other cases also to see the deconfine- 
ment of a boundary-layer giving rise to a shock in the 
bulk (see reference p})- 

With further increase of a.\ , the density variation from 
from x = end, is still the same as before up to part 
(3c) along the /?2+ curve except that now the density ap- 
proaches closer to pj^ along (4c) like path. Finally for a 
given ati, the density reaches the value p*^. After this, 
the boundary condition is satisfied through a depleted 
boundary-layer represented by vertical path (4d). With 
further increase of a\, the density cannot now go around 
p*^ to move to p* + branch due to the constraint from the 
stability property. In that case, the only option for the 
density is to proceed along (4c) but move to p\ + branch 
along a vertical line similar to (5d) before reaching the 
point Pxm' There is now a second high-density upward 
shock in the density profile at larger x (line (5d)) with 
the first shock being a low-density one represented by 
the line, (2b). With the increase of cvi, the (5d) type 
vertical line moves to higher values of Co- The bound- 
ary condition at x = 1 is now satisfied by the rest of 
the density profile where the density decreases continu- 
ously along (5e) path to the minimum p*+ and after that 
decreases through a depleted boundary-layer represented 
by the line (4d). Thus, for example, for this value of evi, 
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FIG. 8: (a) The trajectory of the density on Co — pi plane. 
Five possible trajectories are shown. These trajectories are 
distinguished by different numbers. Alphabetical sequence 
represents the variation of the density along increasing x. (b) 
Plot of the density profiles pi(x) for various a\ with 71 = 
— .28. (c) Zoomed views of the of the boundary-layers of the 
density profiles in (a) near x = 1. 



we see the following parts in the density profile as we 
move along the density profile from its x = end. (i) 
A continuously varying density profile that satisfies the 
boundary condition p(x = 0) = a\. (ii) This is followed 
by a low-density upward shock of maximum height con- 
necting pu = pIq and pi r = 0. (hi) Beyond this shock, 
there is again a continuously increasing part, (iv) This is 
followed by an upward, high-density shock, (v) Beyond 
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this high-density shock, the density decreases continu- 
ously to p*+. (vi) The last part is a particle-depleted 
boundary layer, that saturates to p*{^ for x < 1 and sat- 
isfies the boundary condition at pi(x = 1) = 71. All 
these features can be verified from the density profile in 
figure (Hi). 

If ct\ is increased further, the upward high-density 
shock (vertical lines like (5d)) will move towards higher 
Cq values. For certain ai, the low and the high-density 
shocks merge and there is a big symmetric shock with 
Co = 0. Beyond this a±, the big upward shock still per- 
sists and this is followed by a continuously varying part 
(similar to (lc)) along which the density decreases and 
approaches p^L- The boundary condition is again satis- 
fied by the particle-depleted boundary-layer represented 
by (4d). 



V. GENERAL PREDICTIONS 

Since the boundary-layers or shocks are the special 
features through which the density profiles are distin- 
guished, our general predictions are more on the kind of 
shocks or boundary-layers that can be seen under various 
boundary conditions. 



A. Shocks and boundary-layers 

Downward shock in the bulk or a particle-depleted 
boundary-layer at x = 1: Either of these features appears 
whenever the density profile decreases through a jump 
discontinuity from the P2 + branch to the p\_ branch or 
from to The condition for this is 71 < p*{^ n . 

The value of u\ is somewhat flexible since it is possible 
to see these features both for ol\ positive or negative. 

Upward symmetric shock in the bulk or depleted 
boundary-layer at x = 0: This is seen whenever the den- 
sity profile jumps from p\_ to p 1+ . This happens for var- 
ious combinations of a\ and 71 such as a\ < or otx > 
with 71 > pI^ or 71 < p*+. At x = 0, the density profile 
may start with a continuously varying part followed by 
a symmetric large shock, or it can satisfy the boundary 
condition at x = with the help of a boundary-layer. In 
both cases, the discontinuity in the density corresponds 
to a discontinuous jump from p\_ to p\ + . 

The presence or absence of a boundary-layer at x = 
is specified completely by the value of Co at x — 0. Let us 
assume that p\ = ol\ line intersects the curve on figure Q 
at a value Co = Co(ai). A condition as Co(ai) = Co (a; = 
0), would mean a continuously varying density near x = 
0. If these two values are unequal, it would imply the 
presence of a boundary-layer. For example, for ot\ > 
p*^, a particle- rich or a particle-depleted boundary- layer 
appears if Co(ai) > Cq{x = 0) and Co(ai) < Cq(x = 0), 
respectively. However, it is important to pay attention to 
certain situations which are forbidden due to the stability 



properties. For example, if ot\ < p*~ , a boundary-layer 
with Co (2 = 0) < Co(ai) is not possible. 

Double shock: In this case, the density profile has both 
high and low-density upward shocks with the low-density 
shock having maximum possible height, p*g. In order to 
have a high-density shock, the lower end of the high- 
density shock must be on the P2 + branch. The density 
can reach this branch only via p\ = point. The only 
way the density can reach the p\ = point is through a 
low-density shock represented by the Co = line across 
the negative lobe. A low-density shock representing a 
jump across the negative lobe along Co = line has the 
maximum possible height. 

A density-profile with double shock may appear for 
ai < p\~ n and 71 < p*+ or 71 > pj+. In case of 71 > p*+, 
the density after the high-density shock varies continu- 
ously along p* + branch to satisfy the boundary condition 
at x = 1. For 71 < p*^, the second shock is possible for 
some ot\. In this case, after reaching the p 1+ branch, the 
density decreases till p*^ and then decreases further as a 
depleted boundary-layer at x = 1 to satisfy the boundary 
condition. 

It is interesting to note that although it is possible to 
have a profile with only a low-density shock, the same 
with a single high-density shock is never possible. The 
flow behavior suggests that a high-density shock has to be 
always accompanied by a low-density shock of maximum 
height. 

Boundary-layer at x = 1: As in the case of a boundary- 
layer at x = 1, it is also possible to specify the conditions 
for a boundary-layer at x = 1 by comparing the value of 
Cq(x = 1) with Co(7i). In general, a boundary-layer will 
appear at x = 1 if these two values of Co are different. 
As an example, a downward boundary-layer for < 
71 < appears if C (x = 1) < C (7i). 



B. Saturation of the shock 

From equation (|13|) . we find that near the saturation 
to a bulk density pu>, the slope of the boundary-layer is 
given by 

dSpi J r I u 2 

~W = ( — - 2 P i» )PlbSpi (22) 

where it is assumed that the boundary-layer density is 
6pi away from the saturation value, pi&. This shows, 
that the saturation of the boundary-layer to the bulk 
is in general exponential except for three special points. 
The saturation is of power-law kind, if pi& = or 
Pu = Piti = ^(u) 1 ^ 2, length scale associated with 
the exponential approach of the shock to the bulk den- 
sity diverges as the bulk density approaches these spe- 
cial values. As discussed in subsection III CI the crit- 
ical points, correspond to special boundary conditions 
(a c ,7 c ) at which the shock height across the positive or 
negative lobe reduces to zero. Therefore, the approach to 
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the critical point is associated with the continuous van- 
ishing of the shock height along with the divergence of 
the length scale over which the shock saturates to the 
bulk. 



VI. SUMMARY 

Here, we have considered an asymmetric simple ex- 
clusion process of interacting particles on a finite, one 
dimensional lattice. These particles have mutual repul- 
sion in addition to the exclusion interaction. Apart from 
the hopping dynamics of the particles, the model also has 
particle attachment-detachment processes, which lead to 
particle non-conservation in the bulk. Such processes are 
known to exhibit boundary-induced phase transitions for 
which the tuning parameters are the boundary densities, 
a and 7. In different phases, the average particle den- 
sity distributions across the lattice have distinct shapes 
with various types of discontinuous jumps from one den- 
sity value to another. Here, we carry out a phase-plane 
analysis for the boundary-layer differential equation to 
understand how the fixed-points of the boundary-layer 
equation and their flow properties determine the shape 
of the entire density profile under given boundary con- 
ditions. Such a fixed-point analysis has been extremely 
useful in understanding the phases and phase transitions 
of particle conserving models for which the constant bulk 
density values in different phases are given by the physi- 
cally acceptable fixed-points of the boundary-layer equa- 
tion. In addition, the number of steady-state phases, 
the nature of the phase-transitions, the locations of the 
boundary-layers can be obtained analytically from the 
phase-plane analysis of the boundary-layer equation. The 
present work provides a generalization of the method to 
a particle non-conserving process. 

To apply this method, we have considered the hydrody- 
namic limit of the statistically averaged master equation 
describing the particle dynamics. The hydrodynamic 
equation, describing the time evolution of the average 
particle density, looks like a continuity equation supple- 
mented with the particle non-conserving terms. The cur- 
rent contains the exactly known hopping current and a 
regularizing diffusive current part. The boundary-layer 
equation, which is the main focus of this work, can be 
obtained from the particle conserving part of the hydro- 
dynamic equation. For convenience, we use p± = 2p — 1 
for the boundary-layer equation. p\ is related to the de- 
viation from p —1/2 (half filled case). 

It is found that the fixed-points, p*, of the boundary- 
layer equation are determined in terms of a parameter 
Co related to the excess current measured from p = 1/2 
(half- filled case) . Since the fixed-points are dependent on 
Co , one can plot the physically acceptable fixed-points as 
a function of Co on the Co — p* plane. In the steady- 
state, the constancy of the current across a shock or 
a boundary-layer implies that such objects can be rep- 
resented by a fixed value of Cq. The boundary-layers 



or shocks of the density profiles are represented by the 
constant-Co lines on this Co — p\ plot. The densities 
at which the constant-Co line intersects the fixed-point 
branches are the densities to which the shock or the 
boundary-layer saturates. The discontinuous change of 
the density has to be consistent with the stability prop- 
erties of the fixed-points. For given values of a and 7, 
we can start from the x = end of the density profile 
and find out how the density can change along the pro- 
file as it proceeds to satisfy the boundary condition at 
i = l. This density variation along the density profile 
can be conveniently marked on the Co — p\ plot to see its 
consistency with the flow properties of the fixed-points. 
Our approach does not give any information about the 
location of a shock since it does not involve the details 
of the bulk part of the profile. Instead, it is found that 
the conserved quantity, Co, plays an important role in 
deciding the shape of the density profile. 

The emphasis of our approach is on the boundary-layer 
equation which appears to control the shape of the en- 
tire density profile. Particle non-conserving processes are 
not important for the boundary layers. This simplicity 
allows us not only to classify different kinds of density 
distributions, but also to gain more physical insight as 
why some features of the density profile arc evident un- 
der certain boundary conditions. Some of these features 
are mentioned in the list below. 

(a) When a density profile has two shocks, the low- 
density shock is of maximum possible height. For given 
values of the interaction parameters, the height of the 
low-density shock can be obtained explicitly. 

(b) It is possible to have a low-density shock alone 
in the profile but a high-density shock has to be always 
accompanied by a low-density shock of maximum height. 

(c) A downward shock is produced by the deconfinc- 
ment of a downward boundary-layer at x = 1. The con- 
dition on 7 for seeing a downward shock or a downward 
boundary-layer at x = 1 can be precisely specified. 

(d) The symmetric two peak structure of the current 
as a function of the particle density is responsible for a 
symmetric two lobe structure of the fixed-points drawn 
on Co — p\ plane. The flow behavior of the fixed-points 
around the two lobes are asymmetric. This is the reason 
why the two critical points in the phase diagram are not 
symmetrically related to each other. This asymmetry is 
reflected in the shapes of the density profiles near these 
critical points. 

(e) For a given boundary condition, a density profile 
with only one boundary-layer and no shock can be fully 
specified by the value of Co at this end. 

In addition to these issues, this analysis also provides 
quantitative predictions regarding the heights of different 
kinds of shocks and their approach to the bulk along with 
the length scale associated with it. 
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